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. Abstract. We study the role of viscosity in the early stages of relativistic heavy ion 

I collisions. We investigate the effects of viscosity on the chemical equilibration of a 

parton gas. In the presence of viscosity the lifetime of the system is increased. The 
temperature as well as the parton fugacities evolves more slowly compared to ideal 
fluid dynamics. 



1. Introduction 



>■ ' Ultra-relativistic nucleus-nucleus collisions probe the properties of nuclear matter under 



extreme conditions pQ. Lattice quantum chromodynamics (QCD) calculations [2] predict 



CN ■ that ordinary nuclear matter undergoes a phase transition to quark gluon plasma (QGP). 
CO ' 

An important question is whether the high-energy density matter formed in ultra- 
relativistic nuclear collisions lives sufficiently long enough to reach thermodynamical 
^ ■ equilibrium. That is does the matter reaches thermal, mechanical and chemical 
equilibrium? In this work we assume that the matter reaches thermal and mechanical 
- equilibrium after proper time Tq. We do not, however, assume that the matter is in 

^ ■ chemical equilibrium. Under these assumptions, and given initial values for temperature 



and the quark, antiquark and gluon number densities, we can then employ fluid dynamics 
to study the subsequent evolution of the kinetically equilibrated quark-gluon phase, 
coupled to the rate equations which determine the chemical composition of the system 
far away from equilibrium. This problem has been studied previously in [TTl[T2l[T3] using 
ideal fluid, in [H] using flrst order theory of dissipative fluid dynamics, and recently in 
[15] using parton cascade. 

Ideal (Euler) fluid dynamics has been useful in describing most of the observables 
at RHIC [m [17]. In the early stages of heavy ion collisions, non-equilibrium effects 
play a dominant role. A complete description of the dynamics of heavy ion reactions 
needs to include the effects of dissipation through non-equilibrium/dissipative fluid 
dynamics. As is well-known [HI [10], second order (or extended) theories (which are 
hyperbohc and causal) of dissipative fluids due to Grad [5], Miiller [6j, and Israel and 
Stewart [7] were introduced to remedy some undesirable features such as acausality. 
It seems appropriate therefore to resort to hyperbolic theories instead of flrst order 
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theories (which are parabohc) in describing the dynamics of heavy ion colhsions. First 
order theories are due to Eckart [3] and to Landau and Lifshitz |1] and they lead to 
Navier-Stokes-Fourier (NSF) equations which may be acausal. In addition the first 
order theories do not have well-posed initial value problem. 



2. Viscous hydrodynamics and chemical equilibration 

In this work we extend the work of [8] to include chemically non-equilibrium effects. For 
partially equilibrated plasma of massless particles the equation of state can be written 

as [Tl] 

e = 3p=[a2\g + b2{\ + \g)]T^ , (1) 

where 02 = Svr^/lS, 62 = 7Tc'^Nf/A0 with Nf being the number of quark flavours, and 
the Aj are the parton fugacities defined through 



n, 



9 - K^7^ = ^9 = ^q^f > (2) 

where the n^'^ are the equilibrium parton densities 

n^" = a,T^ = b,T' = nf , (3) 

where oi = (16/7r^)^(3)) and 61 = (9/27r^)^(3). The shear viscosity coefficient is given 
by 

V = \Vq + KVq , (4) 

where the shear viscosity coefficients for the quarks and gluons are given by [20l [2T] 

V, = b,T\ r^g = b,T' , (5) 

where bg = 0.82(a^ ln(l/as))""'^ and bg = 0.20(a^ ln(l/as))^"'^ with being the strong 
coupling constant. We take = 0.4 throughout this analysis, unless otherwise stated. 

In the absence of chemical equilibrium we need the master equations for the 
evolution of parton densities. We consider only the dominant reactions gg ^ ggg, and 
99 QQ- For longitudinal boost invariant longitudinal flow under the assumed equation 
of state and transport coefficients the energy equation and shear pressure [HI [H] become 

Xg + bjXg + X,) t 4 1 ^ 

Xg + b{Xg + q) ^ T^St [a,Xg + b2{Xg + Xg)]T^T ' 

^1 



d£ 2 [a2A, + &2(A, + A,-)] ^^ ^ _ f ^ 

dr 9 [Xgbg + Xgbg] 2 



T 



T Xg + bjXg + X,) 
T ^ A, + 6(A, + q) 



8 [a2Xg + b2{Xg + Xg)]T' 



27 T 

and are coupled to the master equations for the fugacities [HI [121 [13] 
A„ f 1 „ _ _ / A„A 



(7) 



^ + 3- + - = i?3(l-A,)-2i?2(^l-^) , (8) 



Viscous Hydrodynamics 



3 



where b = = 21A^j/64 with Nf being the number of quark flavours The reaction 

rates R2 and R3 are given by [11] 

R2 ~ 0.24Ar/a^AgT ln(1.65/a,Ag) , R3 - 2.1 a^T ^2\g - A2 . (10) 

We consider initial conditions relevant for RHIC: tq = 0.25 fm/c, Tq = 0.66 GeV, 
$0 = po/5, A30 = 0.34, and A^o = 0.064; and for LHC: tq = 0.25 fm/c, Tq = 1.0 GeV, 
"^•0 = Po/5, Ago = 0.43, and Ago = 0.082. Here po is the initial pressure. These sets of 
initial conditions are motivated in [121 E] • 

3. Results 

In Fig. [Tl^a,b) we show the well known results (see ^) of the temperature evolution. 
The ideal fluid dynamics approximation leads to faster cooling. Due to the reduction 
of longitudinal pressure, less work is done to the expansion and hence the slow cooling 
in the presence of viscosity. However the flrst order theory even predicts heating during 
the expansion stage. This is in contradiction to the energy conservation laws. Also the 
flrst order theory will overestimate the freeze-out temperatures. This in turn might lead 
to wrong conclusions about the observables. On the other hand the second order theory 
does not have these undesirable features that are exhibited by flrst order theory. 





1 [fm/c] T [fm/c] 

Figure 1. Time evolution of temperature for (a) RHIC and (b) LHC initial conditions. 
The curves are for ideal fluid (solid), dotted first order theory (dotted) and second order 
theory (dashed). 

In Figs. [2](a,b) and [3](a,b) we show time evolution of the parton fugacities. In 
the presence of viscosity the parton viscosities evolve more slowly towards chemical 
equilibration. This will have considerable effects on the observables such as strangeness 
production. In the early stages of the expansion of the system one sees the difference 
between flrst order theory and second order theory. This difference is investigated in 
more detail [18j . 
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Figure 2. Time evolution of gluon 
fugacity for (a) RHIC and (b) LHC 
initial conditions. 



Figure 3. Time evolution of quark 
fugacity for (a) RHIC and (b) LHC 
initial conditions. 



4. Conclusions 

We have investigated the effects of shear viscosity on the chemical equilibration of the 
parton system in relativistic nuclear collisions. Due to slow cooling of the system in the 
presence of viscosity chemical equilibration is slowed. The effects of transverse expansion 
[9] and of the mass of strange quark in the chemically non-equilibrated viscous system 
are being studied and will be published somewhere |18j . 
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